Polarization Saturation in Multilayered Interfacial Ferroelectrics.

Van der Waals polytypes of broken inversion and mirror symmetries  have been recently shown to exhibit switchable electric polarization even at the ultimate two-layer thin limit. Their out-of-plane polarization has been found to accumulate in a ladder-like fashion with each successive layer, offering 2D building blocks for the bottom-up construction of 3D ferroelectrics. Here, it is demonstrated experimentally that beyond a critical stack thickness, the accumulated polarization in rhombohedral polytypes of molybdenum disulfide saturates. The underlying saturation mechanism, deciphered via density functional theory and self-consistent Poisson-Schrödinger calculations, point to a purely electronic redistribution involving: 1. Polarization-induced bandgap closure that allows for cross-stack charge transfer and the emergence of free surface charge; 2. Reduction of the polarization saturation value, as well as the critical thickness at which it is obtained, by the presence of free carriers. The resilience of polar layered structures to atomic surface reconstruction, which is essentially unavoidable in polar 3D crystals, potentially allows for the design of new devices with mobile surface charges. The findings, which are of general nature, should be accounted for when designing switching and/or conductive devices based on ferroelectric layered materials.

Effect of Interlayer Bonding on Superlubric Sliding of Graphene Contacts: A Machine-Learning Potential Study.

Surface defects and their mutual interactions are anticipated to affect the superlubric sliding of incommensurate layered material interfaces. Atomistic understanding of this phenomenon is limited due to the high computational cost of ab initio simulations and the absence of reliable classical force-fields for molecular dynamics simulations of defected systems. To address this, we present a machine-learning potential (MLP) for bilayer defected graphene, utilizing state-of-the-art graph neural networks trained against many-body dispersion corrected density functional theory calculations under iterative configuration space exploration. The developed MLP is utilized to study the impact of interlayer bonding on the friction of bilayer defected graphene interfaces. While a mild effect on the sliding dynamics of aligned graphene interfaces is observed, the friction coefficients of incommensurate graphene interfaces are found to significantly increase due to interlayer bonding, nearly pushing the system out of the superlubric regime. The methodology utilized herein is of general nature and can be adapted to describe other homogeneous and heterogeneous defected layered material interfaces.

Graphene nanoribbons grown in hBN stacks for high-performance electronics.

Van der Waals encapsulation of two-dimensional materials in hexagonal boron nitride (hBN) stacks is a promising way to create ultrahigh-performance electronic devices1-4. However, contemporary approaches for achieving van der Waals encapsulation, which involve artificial layer stacking using mechanical transfer techniques, are difficult to control, prone to contamination and unscalable. Here we report the transfer-free direct growth of high-quality graphene nanoribbons (GNRs) in hBN stacks. The as-grown embedded GNRs exhibit highly desirable features being ultralong (up to 0.25 mm), ultranarrow (<5 nm) and homochiral with zigzag edges. Our atomistic simulations show that the mechanism underlying the embedded growth involves ultralow GNR friction when sliding between AA'-stacked hBN layers. Using the grown structures, we demonstrate the transfer-free fabrication of embedded GNR field-effect devices that exhibit excellent performance at room temperature with mobilities of up to 4,600 cm2 V-1 s-1 and on-off ratios of up to 106. This paves the way for the bottom-up fabrication of high-performance electronic devices based on embedded layered materials.

Nanotube Slidetronics.

One-dimensional slidetronics is predicted for double-walled boron-nitride nanotubes. Local electrostatic polarization patterns along the body of the nanotube are found to be determined by the nature of the two nanotube walls, their relative configuration, and circumferential faceting modulation during coaxial interwall sliding. By careful choice of chiral indices, chiral polarization patterns can emerge that spiral around the nanotube circumference. The potential usage of the discovered slidetronic effect for low-dimensional nanogenerators is briefly discussed.

Anisotropic Interlayer Force Field for Group-VI Transition Metal Dichalcogenides.

An anisotropic interlayer force field that describes the interlayer interactions in homogeneous and heterogeneous interfaces of group-VI transition metal dichalcogenides (MX2, where M = Mo, W, and X = S, Se) is presented. The force field is benchmarked against density functional theory calculations for bilayer systems within the Heyd-Scuseria-Ernzerhof hybrid density functional approximation, augmented by a nonlocal many-body dispersion treatment of long-range correlation. The parametrization yields good agreement with the reference calculations of binding energy curves and sliding potential energy surfaces. It is found to be transferable to transition metal dichalcogenide (TMD) junctions outside of the training set that contain the same atom types. Calculated bulk moduli agree with most previous dispersion-corrected density functional theory predictions, which underestimate the available experimental values. Calculated phonon spectra of the various junctions under consideration demonstrate the importance of appropriately treating the anisotropic nature of the layered interfaces. Considering our previous parametrization for MoS2, the anisotropic interlayer potential enables accurate and efficient large-scale simulations of the dynamical, tribological, and thermal transport properties of a large set of homogeneous and heterogeneous TMD interfaces.

Ionotronics for reverse actuation.

In the midst of an ongoing energy crisis, the search for new methods of energy harvesting has never been more important. Here we explore, analyse and discuss principles of ionotronic reverse-actuator devices based on the effect of double-layer charging. The designs that we consider in this paper operate based on a common principle - using external mechanical work, which would otherwise be wasted, to produce changes in the contact area of electrode and electrolyte, translated into the time variation of the double-layer capacitance. Periodic variation of capacitance, when connected to a reference voltage source, produces alternating electric current through a load. This concept is not new and in some forms was realised in the early works of Boland, Krupenkin and several papers of our group. The goal of the present paper is to build a comprehensive analytical platform for a description of operation of such devices in terms of materials, generated power as a function of the frequency of variation of applied force, electrical load, and other factors; the understanding of which allows us to optimise these systems and navigate their construction. The first design, discussed in the paper, is based on flat electrodes. It is the simplest one and, as such, helps elucidate some key factors determining power generation. While being easy to realise experimentally, it generates relatively low power, even when optimised. The second design, based on microporous electrodes is more sophisticated and allows a much larger power harvest. The results are also compared to the recently proposed capacitive rotor device. The developed theory is set to capture the key factors that determine the functioning of the considered reverse actuators. The structures under study are matched to fit into the sole of a shoe and produce power from walking and running. However, they can also be scaled-up to larger operating systems and various external loads.

Electric-field frictional effects in confined zwitterionic molecules.

We theoretically explore the effect of a transverse electric field on the frictional response of a bi-layer of packed zwitterionic molecules. The dipole-moment reorientation promoted by the electric field can lead to either stick-slip or smooth sliding dynamics, with average shear stress values varying over a wide range. A structure-property relation is revealed by investigating the array of molecules and their mutual orientation and interlocking. Moreover, the thermal friction enhancement previously observed in these molecules is shown to be suppressed by the electric field, recovering the expected thermolubricity at large-enough fields. The same holds for other basic tribological quantities, such as the external load, which can influence friction in opposite ways depending on the strength of the applied electric field. Our findings open a route for the reversible control of friction forces via electric polarization of the sliding surface.

Orientational Ordering in Nano-confined Polar Liquids.

Water and other polar liquids exhibit nanoscale structuring near charged interfaces. When a polar liquid is confined between two charged surfaces, the interfacial solvent layers begin to overlap, resulting in solvation forces. Here, we perform molecular dynamics simulations of polar liquids with different dielectric constants and molecular shapes and sizes confined between charged surfaces, demonstrating strong orientational ordering in the nanoconfined liquids. To rationalize the observed structures, we apply a coarse-grained continuum theory that captures the orientational ordering and solvation forces of those liquids. Our findings reveal the subtle behavior of different nanoconfined polar liquids and establish a simple law for the decay length of the interfacial orientations of the solvents, which depends on their molecular size and polarity. These insights shed light on the nature of solvation forces, which are important in colloid and membrane science, scanning probe microscopy, and nano-electrochemistry.

Theory and Simulations of Ionic Liquids in Nanoconfinement.

Room-temperature ionic liquids (RTILs) have exciting properties such as nonvolatility, large electrochemical windows, and remarkable variety, drawing much interest in energy storage, gating, electrocatalysis, tunable lubrication, and other applications. Confined RTILs appear in various situations, for instance, in pores of nanostructured electrodes of supercapacitors and batteries, as such electrodes increase the contact area with RTILs and enhance the total capacitance and stored energy, between crossed cylinders in surface force balance experiments, between a tip and a sample in atomic force microscopy, and between sliding surfaces in tribology experiments, where RTILs act as lubricants. The properties and functioning of RTILs in confinement, especially nanoconfinement, result in fascinating structural and dynamic phenomena, including layering, overscreening and crowding, nanoscale capillary freezing, quantized and electrotunable friction, and superionic state. This review offers a comprehensive analysis of the fundamental physical phenomena controlling the properties of such systems and the current state-of-the-art theoretical and simulation approaches developed for their description. We discuss these approaches sequentially by increasing atomistic complexity, paying particular attention to new physical phenomena emerging in nanoscale confinement. This review covers theoretical models, most of which are based on mapping the problems on pertinent statistical mechanics models with exact analytical solutions, allowing systematic analysis and new physical insights to develop more easily. We also describe a classical density functional theory, which offers a reliable and computationally inexpensive tool to account for some microscopic details and correlations that simplified models often fail to consider. Molecular simulations play a vital role in studying confined ionic liquids, enabling deep microscopic insights otherwise unavailable to researchers. We describe the basics of various simulation approaches and discuss their challenges and applicability to specific problems, focusing on RTIL structure in cylindrical and slit confinement and how it relates to friction and capacitive and dynamic properties of confined ions.

Interlayer Registry Dictates Interfacial 2D Material Ferroelectricity.

We discover that the complex ferroelectric response of layered materials toward interlayer sliding is fully dictated by the interlayer lattice registry. Importantly, the entire sliding polarization landscape of two-dimensional (2D) layered material interfaces is fully described via a simple and intuitive geometric measure, termed the polarization registry index (PRI), that quantifies the degree of interlayer commensurability. Beyond the understanding of the fundamental origin of 2D ferroelectricity, the developed tool also provides highly efficient characterization and rationalization of existing experimental and computational evidence of 2D interfacial ferroelectricity, as well as the prediction of emergent controllable polarization in new noncentrosymmetric layered systems.

Velocity Dependence of Moiré Friction.

Friction force microscopy experiments on moiré superstructures of graphene-coated platinum surfaces demonstrate that in addition to atomic stick-slip dynamics, a new dominant energy dissipation route emerges. The underlying mechanism, revealed by atomistic molecular dynamics simulations, is related to moiré ridge elastic deformations and subsequent relaxation due to the action of the pushing tip. The measured frictional velocity dependence displays two distinct regimes: (i) at low velocities, the friction force is small and nearly constant; and (ii) above some threshold, friction increases logarithmically with velocity. The threshold velocity, separating the two frictional regimes, decreases with increasing normal load and moiré superstructure period. Based on the measurements and simulation results, a phenomenological model is derived, allowing us to calculate friction under a wide range of room temperature experimental conditions (sliding velocities of 1-104 nm/s and a broad range of normal loads) and providing excellent agreement with experimental observations.

Cumulative polarization in conductive interfacial ferroelectrics.

Ferroelectricity in atomically thin bilayer structures has been recently predicted1 and measured2-4 in two-dimensional materials with hexagonal non-centrosymmetric unit-cells. The crystal symmetry translates lateral shifts between parallel two-dimensional layers to sign changes in their out-of-plane electric polarization, a mechanism termed 'slide-tronics'4. These observations have been restricted to switching between only two polarization states under low charge carrier densities5-12, limiting the practical application of the revealed phenomena13. To overcome these issues, one should explore the nature of polarization in multi-layered van der Waals stacks, how it is governed by intra- and interlayer charge redistribution and to what extent it survives the addition of mobile charge carriers14. To explore these questions, we conduct surface potential measurements of parallel WSe2 and MoS2 multi-layers with aligned and anti-aligned configurations of the polar interfaces. We find evenly spaced, nearly decoupled potential steps, indicating highly confined interfacial electric fields that provide a means to design multi-state 'ladder-ferroelectrics'. Furthermore, we find that the internal polarization remains notable on electrostatic doping of mobile charge carrier densities as high as 1013 cm-2, with substantial in-plane conductivity. Using density functional theory calculations, we trace the extra charge redistribution in real and momentum spaces and identify an eventual doping-induced depolarization mechanism.

Electrotunable friction with ionic liquid lubricants.

Room-temperature ionic liquids and their mixtures with organic solvents as lubricants open a route to control lubricity at the nanoscale via electrical polarization of the sliding surfaces. Electronanotribology is an emerging field that has a potential to realize in situ control of friction-that is, turning the friction on and off on demand. However, fulfilling its promise needs more research. Here we provide an overview of this emerging research area, from its birth to the current state, reviewing the main achievements in non-equilibrium molecular dynamics simulations and experiments using atomic force microscopes and surface force apparatus. We also present a discussion of the challenges that need to be solved for future applications of electrotunable friction.

Catalytic Growth of Ultralong Graphene Nanoribbons on Insulating Substrates.

Graphene nanoribbons (GNRs) with widths of a few nanometers are promising candidates for future nanoelectronic applications due to their structurally tunable bandgaps, ultrahigh carrier mobilities, and exceptional stability. However, the direct growth of micrometer-long GNRs on insulating substrates, which is essential for the fabrication of nanoelectronic devices, remains an immense challenge. Here, the epitaxial growth of GNRs on an insulating hexagonal boron nitride (h-BN) substrate through nanoparticle-catalyzed chemical vapor deposition is reported. Ultranarrow GNRs with lengths of up to 10 µm are synthesized. Remarkably, the as-grown GNRs are crystallographically aligned with the h-BN substrate, forming 1D moiré superlattices. Scanning tunneling microscopy reveals an average width of 2 nm and a typical bandgap of ≈1 eV for similar GNRs grown on conducting graphite substrates. Fully atomistic computational simulations support the experimental results and reveal a competition between the formation of GNRs and carbon nanotubes during the nucleation stage, and van der Waals sliding of the GNRs on the h-BN substrate throughout the growth stage. This study provides a scalable, single-step method for growing micrometer-long narrow GNRs on insulating substrates, thus opening a route to explore the performance of high-quality GNR devices and the fundamental physics of 1D moiré superlattices.

Interlayer Registry Index of Layered Transition Metal Dichalcogenides.

Inspired by the fascinating electronic properties of twisted transition metal dichalcogenides, we extend the registry index approach to quantify the interlayer commensurability of homogeneous and heterogeneous interfaces of MoS2, WS2, MoSe2, and WSe2. The developed geometric measure provides quantitative information about their sliding energy landscape with vast mechanical and tribological implications. Furthermore, the registry index is highly suitable for characterizing surface reconstruction in twisted transition metal dichalcogenide interfaces that dictates their intricate electronic and ferroelectric properties. The simple and intuitive nature of the registry index marks it as a powerful computational tool for studying the fascinating physical phenomena demonstrated by these materials.

Structural Forces in Ionic Liquids: The Role of Ionic Size Asymmetry.

Ionic liquids (ILs) are charged fluids composed of anions and cations of different size and shape. The ordering of charge and density in ILs confined between charged interfaces underlies numerous applications of IL electrolytes. Here, we analyze the screening behavior and the resulting structural forces of a representative IL confined between two charge-varied plates. Using both molecular dynamics simulations and a continuum theory, we contrast the screening features of a more-realistic asymmetric system and a less-realistic symmetric one. The ionic size asymmetry plays a nontrivial role in charge screening, affecting both the ionic density profiles and the disjoining pressure distance dependence. Ionic systems with size asymmetry are stronger coupled systems, and this manifests itself both in their response to the electrode polarization and spontaneous structure formation at the interface. Analytical expressions for decay lengths of the disjoining pressure are obtained in agreement with the pressure profiles computed from molecular dynamics simulations.

Thermal Friction Enhancement in Zwitterionic Monolayers.

We introduce a model for zwitterionic monolayers and investigate its tribological response to changes in applied load, sliding velocity, and temperature by means of molecular-dynamics simulations. The proposed model exhibits different regimes of motion depending on temperature and sliding velocity. We find a remarkable increase of friction with temperature, which we attribute to the formation and rupture of transient bonds between individual molecules of opposite sliding layers, triggered by the out-of-plane thermal fluctuations of the molecules' orientations. To highlight the effect of the molecular charges, we compare these results with analogous simulations for the charge-free system. These findings are expected to be relevant to nanoscale rheology and tribology experiments of locally-charged lubricated systems such as, e.g., experiments performed on zwitterionic monolayers, phospholipid micelles, or confined polymeric brushes in a surface force apparatus.

Negative tension controls stability and structure of intermediate filament networks.

Networks, whose junctions are free to move along the edges, such as two-dimensional soap froths and membrane tubular networks of endoplasmic reticulum are intrinsically unstable. This instability is a result of a positive tension applied to the network elements. A paradigm of networks exhibiting stable polygonal configurations in spite of the junction mobility, are networks formed by bundles of Keratin Intermediate Filaments (KIFs) in live cells. A unique feature of KIF networks is a, hypothetically, negative tension generated in the network bundles due to an exchange of material between the network and an effective reservoir of unbundled filaments. Here we analyze the structure and stability of two-dimensional networks with mobile three-way junctions subject to negative tension. First, we analytically examine a simplified case of hexagonal networks with symmetric junctions and demonstrate that, indeed, a negative tension is mandatory for the network stability. Another factor contributing to the network stability is the junction elastic resistance to deviations from the symmetric state. We derive an equation for the optimal density of such networks resulting from an interplay between the tension and the junction energy. We describe a configurational degeneration of the optimal energy state of the network. Further, we analyze by numerical simulations the energy of randomly generated networks with, generally, asymmetric junctions, and demonstrate that the global minimum of the network energy corresponds to the irregular configurations.